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Abstract 

We perform numerical simulations of the two and three-dimensional 
spin systems with competing interaction. They describe the model of ran- 
dom surfaces with linear-gonihedric action. The degeneracy of the vacuum 
state of this spin system is equal to d ■ 2 for the lattice of the size 
N d . We observe the second order phase transition of the three-dimensional 
system, at temperature [3 C ~ 0.439322 which almost coincides with (3 C of 
the 2D Ising model. This confirms the earlier analytical result for the case 
when self-interaction coupling constant k is equal to zero. We suggest the 
full set of order parameters which characterize the structure of the vacuum 
states and of the phase transition. 



U.Crete TH 18/95 
February 1995 



1. 

It is well known that the partition function of three-dimensional Ising ferro- 
magnet is equivalent to the summation over closed surfaces with an action which 
is proportional to the surface area H |], || . In two and three dimensions the Ising 
model undergoes a second order phase transition and at the critical point, in 
two- dimensions, corresponds to free fermions and probably to interacting string 
theory in three dimensions [|, [7], §, §, |l(|. 

The aim of this article is to investigate the phase structure of the two- and 
three-dimensional spin systems introduced in 0. In three dimensions the par- 
tition function of this spin system can be represented as a sum over closed self- 
intersecting random surfaces with an action which is proportional to the linear 



size of the surfaces [11, 12, 13 



We observe the second order phase transition of the system with linear- 
gonihedric action in three dimensions. This numerical result encourages the hope 
that the system may describe the propagation of an almost free fermionic string 
||. This hope is based on the fact that the the 3D gonihedric model can be 
rewritten as a model of propagating loops with the interaction term which is 
proportional to the overlapping length, while 3D Ising model can be rewritten in 
the same way, but with the interaction term which is proportional to overlapping 
area, and therefore in the latter case the interaction is much stronger 0. 

In the next sections we describe the two- and three-dimensional models with 
gonihedric action. We stress that the vacuum is strongly degenerated and that 
the rate of degeneracy depends on the self-intersection coupling constant k. If 
k 7^ 0, the degeneracy of the vacuum state is equal to d ■ 2 N for the lattice 
of the size N d and is equal to 2 dN when k = 0. The last case is a sort of 



super symmetric point in the space of gonihedric Hamiltonians [14], The rate 
of degeneracy is between Ising model Z 2 and Wegner gauge spin systems ■ 
In two dimensions the system is in completely disordered regime because of 
the tunneling phenomenon between the vacuum states through the instantons. 
The instantons remove the degeneracy of the ground state and the symmetry 
restoreds. 

In three dimensions the vacua are well separated and we have different sym- 
metry breaking phases. We observed the second order phase transition at j3 c ~ 
0.439322. We suggest the full set of order parameters which characterize the 
structure of the vacuum states and the phase transition. At the critical point the 
correlation length became infinite and the theory may have a continuum limit. 

2. 
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The Hamiltonian of the system in two dimensions has the form jlL IT3 



k 1 — k 

Hgonimetric = ~ k ^f+S + ^ ^ Yl ^f+S '?+S+0° r+/3 W 

The low temperature expansion of the partition function 

Z(P) = £ exp (-/3 H 2 g d onimetric ) (2) 
W 

can be represented as a sum over random walks or paths with the energy of the 
path which is proportional to the number of right corners of the path. In the 
self-intersection points the interaction energy is equal to 4k and therefore the 
constant k is called an intersection coupling constant jl3[ [14]] • The energy of a 
given path P is equal to 



E path = k{P) = right(P) + Ak ■ inter(P) (3) 

where right(P) is the number of right angles of the path and inter(P) is the 
number of its self-intersection points. E pat h is a scale invariant quantity: if we 
enlarge the contour let's say A times, then E path does not change. 
Comparing this system with 2d Ising ferromagnet 

H I sing = ~ r+d ( 4 ) 

f,a 

one can see that the energy of paths in the low-temperature expansions of the 
partition function is proportional to the length of the path 

E pa t h = 1{P) = Unks{P) (5) 

where links(P) is the number of the links of the closed path. 

The intersection coupling constant k defines the intensity of self-intersections 
of the paths and at the same time defines the rate of degeneracy of the vacuum 
states. When k ^ 0, then the vacuum is degenerated by the factor d ■ 2 N for 
the lattice of the size N d , where d is the dimension of the lattice. The degeneracy 
increases when k = and is equal to 2 dN . This means that the point k = in the 
"space" of Hamiltonians corresponds to a system with a higher symmetry. The 
rate of degeneracy of the vacuum is between the Ising ferromagnet and Gauge-spin 
models which is Z 2 and Z^ d correspondingly. 
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In the case when the intersection coupling constant is equal to zero (k = 0) 
the Hamiltonian (1) and the partition function reduce to the form 

H gonimetric = ~ 9 a r®r+aO ~f?+a+/3° ~r+0 (P) 
r,a,/3 

Z(0) = £ exp (-(3 Ef onmetnc ) = £ e-" k(P \ (?) 
M {P} 

where k(P) is the total curvature of the path (3). In this particular case k = 

the model can be solved exactly 0. Indeed the high temperature expansion of 

the partition function is: 

Z((3) = £ exp (-(3 Ef ommetnc ) = ( 2 ch(5) N \l + 0(N) + . . .) (8) 
M 

and the free energy is equal to 

-Pf = ln(2ch(3) (9) 

The system is in a completely disordered regime like the one-dimensional Ising 
model with energy density and the specific heat equal to: 

The reason why the system is in the disordered regime is connected with the fact 
that the potential barrier between these vacua is finite like in the ID Ising model 
and the tunneling phenomenon through the instantons removes the degeneracy 
of the ground state and the symmetry restores. 

The question to which we would like to address here is the nature of the phase 
transition in the case k ^ 0. The model with k ^ can be formulated as an eight 



vertex model with the weights [14 



U>1 — 1, U)2 — LU Ak , Co> 3 = U4 = 1 , CJ5 = = UJj = CJg = UJ (11) 

where uo = exp(—(3), and the partition function can be rewritten as a fermionic 
partition function with quadratic and quartic field operators. Perturbation ex- 
pansion around quadratic operators {k = 0) shows that the system is in the 
disordered regime for small values of intersection coupling constant k. 

To confirm this universal behavior of the system for different values of the 
self-intersection coupling constant k one can use the Monte-Carlo simulation of 
the spin system (1). Let us consider the case k = 1 0. 
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n, 



gonimetric 



2d 




r,a,0 



(12) 



We performed a simulation of the system on a 100 2 lattice with periodic boundary 
conditions. We measured the magnetization M. =< a >p, the spin-spin corre- 
lation function and the average energy density. We saw that in two dimensions 
there is no phase transition for any positive temperature in accordance with the 
arguments above. 

In order to measure magnetization we start from a high value of (5 with a 
completely magnetized state M. = 1. While we heat up the system, magnetization 
drops and the remaining fluctuations are around M. = 0. As we cool down 
again, we see that the system does become magnetized spontaneously. This is 
the first similarity encountered between the 2D gonimetric model and ID Ising. 
The second similarity comes from the measurements of the spin-spin correlation 
function G(r) =< <r(0),o-(r) > which was seen to drop very fast with distance 
in the large temperature interval. Our measurements of the critical exponent 
v for the correlation length £ at f3 c = oo indicate that v is close to one. The 
third similarity observed is that typical configurations for high f3 have a simple 
structure: the islands formed by spins are rectangular blocks, the system becomes 
more and more ordered as (3 becomes higher. At high temperature the islands 
have curly boundaries, but still the intersections are very disfavoured. This is a 
consequence of the high energy cost associated with them. 

Again one can explain this disordered phase by analyzing the vacuum struc- 
ture. As we already mentioned in the case of Hamiltonian (1) with k ^ the 
vacuum is 2 • 2 N times degenerated, because flat layers of spins with opposite 
directions have zero interface energy. The vacuum state of the 2D Ising model is 
twice degenerated and the potential barrier between those vacua is infinite, while 
in the gonimetric model (1) the barriers between the vacuum states are finite and 
are equal to four units of energy. This is the reason why the system is always 
in a disordered regime. The kinklike excitations which connect those vacua are 
shown on Fig. 1. 

In summary for the two-dimensional system with gonimetric action we have 
the following properties: (i) The system has strong degeneracy of the vacuum 
states, similar to gauge or spin glass systems, which are separated by finite po- 
tential barriers. There is a tunneling phenomenon between those vacua through 
instantons. (ii) Consequently, the instantons remove the degeneracy of the ground 
state and the symmetry restoreds. The system is in the disordered regime for all 
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non-zero temperatures and for large values of the intersection coupling constant 
k. (iii) The typical configurations of interface are such that the paths of interface 
are self-avoiding. 
3. 

In three dimensions the corresponding Hamiltonian is [[1], [H}] 



(13) 



and the low temperature expansion of the partition function 

Z(f3)=Y.exp(-(3H 3 g d omhednc ) (14) 
M 

can be expressed as a sum over surfaces of interface with the amplitude which is 



proportional to linear size of the surface [11, [L2L O 



Zifl) = E e ~ M{M) > A(M) = £ \j ■ |tt - a h3 \ (15) 

{A/} <i,j> 

where Ajj is the length of the surface edges (on the lattice A is equal to lattice 
size a), and aij is the dihedral angle equal to , n/2 or it . 

The constant k again plays the role of intersection coupling constant and 
defines the rate of degeneracy of the vacuum. If we take periodic boundary 
conditions and k = 0, then the vacuum is degenerated by the factor 3 • 2 N and 
consists of spin layers of different width in the x, y, and z direction. In the same 
way as in two dimensions the system simplifies in the symmetric point where the 
self-intersection coupling constant is zero k = ]T3|, 0] 



H gonihedric — 2 ^ CX r a r+aC r f +3+ 0O r~ + (16) 
f,a,beia 

The vacuum is strongly degenerate 2 N ■ 2 N ■ 2 N times , because now all flat 
layers in all directions form different vacuum states. The system (17) is highly 
symmetric and even allows to construct the dual Hamiltonian (the point k = 
in the space of Hamiltonians) . 

Little is known about phase transition point of the three-dimensional system 
(14) even in supersymmetric point k = (17) . Curvature representation of the 
linear action (16) |15| allows to find equivalent representations of the partition 
function for the system (17) in terms of propagation of the polygon loops or 
strings P in a given direction with the transition amplitude 

exp (-k(P) - 2 ■ l(P)) (17) 
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where l(P) is the length and k(P) is the absolute curvature of P. The interaction 
is proportional to the overlapping length of the loops 

A mt = l(P 1 DP 2 ). (18) 

In the first approximation ignoring the interaction term one can solve the model 
and see that the system undergoes a second order phase transition at (3 C which 
is similar to 2D Ising ferromagnet and that it describes the propagation of an 
almost free string of 2D Ising fermions 0. 

One can understand the nature of this second order phase transition analyzing 
the structure of the vacuum. Indeed the potential barrier between those vacua 
is infinite like in 2D Ising model and is of order N for the lattice of the size iV 3 . 
Therefore one can have well separated phases of the system at low temperature, 
but rigorous prove of there existence is still missing. The main problem is to 
count properly the entropy of the given configurations. 

This symmetry breaking phenomenon cannot be uniquely characterized by an 
order parameter like magnetization, because the number of vacuum states with 
zero magnetization is equal to 

N\ 2 N 

WW " 7n (19) 

and the probability to have zero magnetization dominates. The degenaracy is 
less for the vacuum states with nonzero magnetization. The magnetization does 
not uniquely characterize the system and only partly play the role of an order 
parameter of the system. One can introduce the full set of order parameters to 
characterize this phases. Let as denote the vacuum spin comfigurations of the 
system by a^ivac) where \i = 1, 2, 2 m , then the generalized magnetization 
M M is equal to 

M M =< Y, a?(vac) ■ o ? > (20) 

f 

and uniquely characterize the vacuum statets and the phase transition. The 
internal energy, which is equal to the number of corners can also be used as an 
order parameter. 

It is interesting therefore to have Monte-Carlo simulations of the system in 
three dimensions and to study the nature and the dependence of the phase tran- 
sition as a function of the self-intersection coupling constant k in (14). We will 
consider here the case k = 1 ]T[, when 

Hgonihedric = (J r (J r+a + ^ E {T P <7 f+S+0 (^1) 
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and the spin system includes only the competing ferromagnetic and antiferro- 
magnetic interaction. 

We consider the 3D lattice of size iV 3 = 24 3 and 32 3 . Measuring the energy 
density 

«w> = = W3 wm (22) 

we observe the sharp behavior at temperature (3 C ~ 0.439322, see Fig. 2. The two 
point correlation function drastically increase at critical temperature, see Tab.l. 
As we approach the transition point from higher temperatures the correlation 
function blows- up uniformly in all three directions. Afther critical point, at low 
temperatures, the correlation function is unisotropic and we observe symmetry 
breaking phenomena. All energy histograms obtained near f3 c have only one pick, 
see Fig. 3. The typical spin configurations, while we set near critical point , 
describe large surfaces without self-intersections. 

We came to the conclusion that the system in three dimensions undergoes 
a second order phase transition and therefore may have well defined continuum 
limit. For a rigorous proof of this picture one should have a more detailed study 
of the phenomena near the critical point. We still dont know the critical indices 
of the model to distinguish them from the ones in 3D Ising model. 
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Figure caption 

Fig. 1. (a) one of the vacuum states, (b) excitation - the kink of energy four, 
) new vacuum after transition . 

Fig. 2. Total number of corners as a function of (3 in 32 3 3d lattice, at 
= 1. 

Tab. 1 Corralation function at (5 = 0.439320. 
Fig. 3. Energy histogram at (5 = 0.439322 



10 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9504054vl 
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